Detecting paired and counterflow superfluidity via dipole oscillations 
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We suggest an experimentally feasible procedure to observe counterflow and paired superfluidity 
in ultra-cold atom systems. We study the time evolution of one-dimensional mixtures of bosonic 
atoms in an optical lattice following an abrupt displacement of an additional weak confining poten- 
tial. We find that the dynamic responses of the paired superfluid phase for attractive inter-species 
interactions and the counterflow superfluid phase for repulsive interactions are qualitatively distinct 
and reflect the quasi long-range order that characterizes these phases. These findings suggest a clear 
experimental procedure for their detection, and give an intuitive insight into their dynamics. 
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Rapid progress in experiments with ultra-cold atomic 
mixtures has enabled the study of the rich quantum 
many-body phenomena of strongly correlated multi- 
component systems in a controllable environment A 
central current objective of the ultra-cold atom commu- 
nity is the realization and study of magnetic order, and, 
closely related, of 'J^/[/ '-driven physics. J refers to the 
tunneling energy, and U refers to the interaction strength 
in a Hubbard model, discussed below. The most promis- 
ing choice is to use bosonic atoms, due to the higher de- 
gree of degeneracy that can be achieved in such systems. 
In [3] we identify the regime in which bosonic mixtures in 
optical lattices sustain counterflow superfluidity (CFSF). 
Requiring only repulsive contact interactions, this order 
is characterized by particles binding with holes of the 
other species. Such particle-hole pairing, or anti-pairing, 
leads to non-dissipative counter-flow of the two species 
while the net flow is zero Besides CFSF, the ground 
state can also display paired superfluidity (PSF), for at- 
tractive interactions [3|-[9|- 

We propose to realize CFSF. A key question is what 
maximal value of J^/J7 can be used. CFSF competes en- 
ergetically with the single-particle superfluid (SF) phase. 
In [3] we demonstrate that with U12/U w 0.6... 0.8, 
CFSF is sustained up to J/U w 0.15. U12 is the interac- 
tion strength between the two species, see below. When 
using, say, a mixture of two hyperfine states of ^^Rb in 
an optical lattice this can be achieved by dislocating the 
two species slightly with a magnetic field gradient (lO| . 

In this paper we propose to use dissipationless coun- 
terflow as the experimental signature of CFSF. In ex- 
periments of ultra-cold atoms in optical lattices, trans- 
port properties are often studied by suddenly displacing 
a confining harmonic potential and inducing dipole oscil- 
lations. This type of experiment has been carried out to 
study superfluidity of ID (TTI.[T2| and 3D [l^ Bose gases. 
Here, we report the first many-body simulation to ob- 
serve CFSF 12]. We consider SF, PSF and CFSF states in 
a ID Bose mixture in an optical lattice confined in a weak 
trapping potential [3!,[5l]. Using the quasi-exact numerical 
method of time-evolving block decimation (TEBD) [l3], 
we study the transport properties through the dipole os- 
cillations induced by either a brief or a constant displace- 
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Figure 1: (Color online) Dipole oscillations in the CFSF state 
after a brief displacement of species f, where J/U — 1/14, 
U12/U = 0.6, n/U ^ 10"^ and time t is in units of h/J. 
(a), (b): density distributions of species 1 and 2, respectively. 
At t = 0, each species has a plateau at half-filling near the 
trap center. Then, the trap of species 1 is perturbed by a 
brief displacement as shown in Fig. [2ja). Species 1 and 2 
are seen to begin moving in opposite directions, whereas the 
total density is independent of time as shown in panel (c). (d): 
local current of species 1, /ij , defined in Eq. ([2]). A positive 
current corresponds to motion to the right. The diamond 
pattern in panel (d) refiects the constant speed of fiow. 



ment of the harmonic trapping potential. Finding the 
qualitative features described here experimentally, would 
demonstrate the existence of CFSF order. We note that 
the same physical effects can also be achieved by apply- 
ing magnetic field gradients, and that similar results can 
be expected for higher dimensions. 

We consider a ID two-component Bose Hubbard model 
with a harmonic trap centered at Ca(t), 

^ = E E {-Jiblj + lba.j + h.C.) + r! [j - Ca{t)f n,. 



a=l,2 j 

u 



+ lT"aj(naj - 1) ^ + Ui2^ni 



(1) 



We denote the atomic species with index a, the lattice 
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Figure 2: (Color online) Center of mass oscillations of the two 
species for SF, PSF, and CFSF states induced by a brief trap 
displacement of species 1. (a): trap displacement vs. time, 
(b): evolution of the SF state with parameters J/U — 1/8, 
U12/U = -0.1, n/U = 2.5 X 10"*. The full (dashed) lines 
correspond to species 1 (2). (c): PSF state for J/U = 1/8, 
U12/U = -0.7, n/U ^ 2.5 X 10"^. The species oscillate in 
phase as a result of their pairing order, (d): CFSF state for 
J/U = 1/8, U12/U = 0.7 and Q/U = 10"^ The species 
oscillate out of phase due to the anti-pairing order. 



site with index j, the number of lattice sites TV, and we 
impose hard-wall boundary conditions. In this paper, 
N — 80. The species have the same average filling factor, 
1/ = M/N < 1, where M is the number of particles for 
each species. For all cases reported here, M = 20. We 
also assume that the repulsive intra-species interaction 
U > 0, hopping parameter J > and the spring constant 
f2 are the same for both species. The inter-species inter- 
action is given by Ui2- The operators b^^ j and ba.j are 
the creation and annihilation operators for atoms of type 
a on site j and Uaj — b\ jba,j is the atom number op- 
erator. We assume that the trap centers, Ca{t). are time 
dependent and independently controllable. Their initial 
values are Ca{t = 0) = 40.5. For convenience, we define 
the displacement Da{t) = Ca{t) — Ca(0). All distances are 
in units of the lattice constant d. 

The initial state is the ground state of _ff at t = 
obtained by the TEBD method with imaginary time 
propagation. The time evolution is obtained with real- 
time TEBD propagation with the time step 5t equal to 
dtJ/h = 0.05, where H is the reduced Planck constant. 
Previous works [isl . [l6| have applied TEBD to simulate 
dipole oscillations in ID single-species systems. The sim- 
ulations here use similar parameters as in Ref.[15]. We 
analyze the time evolution of the system by studying the 
spatial density distribution {naj{t)), where the expecta- 
tion value is over the state of the system, and also the 
local fiux or current 



fa,j = l{Jiml, + lK3 - &I,&aj + l), (2) 
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Figure 3: (Color online) Dipole oscillation of the total and 
relative center of mass; initial displacement as in Fig. [2] (a): 
xtot{t) for U/J = 2, 6, 10, 12; only in the cases of U/J = 2 
and 6 do we see significant induced oscillation, (b): Xrei{t) 
for the same U/J. A phase transition from SF to CFSF is 
induced by changing U/J while keeping U12/U and i}/U fixed 
at U12/U = 0.6 and Q/U = 10~^. The transition occurs for 
(7/ J ~ 8. The transition is reflected in the sudden change in 
the behavior of the center of mass motion. Panel (c) shows the 
ratio {hLf/{dJTf as a function of C/12/ J in the CFSF region 
for fixed U12/U, where T is the oscillation period and L is the 
length of the CFSF plateau extracted from the simulations. 
The markers are data from our simulations. The full line is a 
linear fit and the dashed line corresponds to Eq. 



and the center of mass for each species 
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We also define the total (relative) center of mass, Xt{t) = 
[xi{t) + X2{t)] /2 and Xr{t) = xi{t) — x^it). 

We first consider the time evolution after a brief dis- 
placement of the harmonic confinement of species 1. 
Figure [1] shows a compelling example of the counter- 
flow property of the CFSF phase. Particle-hole pair- 
ing in the CFSF state requires that the particle den- 
sity of one species equals the hole density of the other: 
(rij^i) = (1 — rtj,2}- For equal fillings, this means that 
the CFSF phase only occurs at half-filling, u =1/2^. 
In a trapped system, the density of the CFSF phase 
shows a plateau at half-filling near the center of the trap. 
Near the edges, the system is SF. The impulse applied 
to species 1 is generated by the brief displacement shown 
in Fig. IHa). It causes the species to move in opposite 
directions, as shown in Fig.[IJa) and (b). The density of 
species a at a given position oscillates about Ua — 0.5, 
and the density oscillations of the two species differ in 
phase by tt. The flux is reflected at the steep SF edge of 
the atomic cloud. The total density distribution stays un- 
changed and equals one (see Fig.[ljc)). Figure[ljd) shows 
the local current of species 1 as a function of time. At any 
time, the flow occurs in only one direction. Moreover, the 
color image shows characteristic diamond shapes. Within 
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Figure 4: Dipole oscillation of the center of mass of the two 
species after a sudden displacement, (a): trap displacement 
vs. time, (b), (c), (d): time evolution for the SF, PSF and 
CFSF states. The system parameters are the same as in Fig. [5] 
and the full (dashed) lines correspond to species 1 (2). The 
depicted time scale differs between panels. 
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Figure 5: Dipole oscillation of the center of mass after dis- 
placing the two traps in opposite directions, (a): trap dis- 
placement vs. time, (b): center of mass motion for the SF, 
PSF and CFSF states. The full (dashed) line corresponds to 
species 1 (2). (c)-(e): density distributions at tJ/h = 55, 100 
and 400, respectively, for the CFSF state. Parameters used 
for the three states are the same as in Fig. [2] 



each diamond, the current is nearly independent of lat- 
tice site and time, indicating a constant velocity flow. 
This constant flow can be understood from the constant 
total density of the CFSF state. As the energy due to 
the trapping potential at each site only depends upon the 
total density there, there is no potential energy cost in 
changing the local relative density and the observed flow 
moves freely within the CFSF plateau. 

In Fig. [21 we show an example of the center of mass 
motion, Xait), for SF, PSF, and CFSF initial states after 
applying a small impulse to species 1. For the SF state. 
Fig. mb) shows that only species 1 is excited immedi- 
ately after the impulse. Oscillatory motion of species 2 is 
only induced gradually as a result of the weak attraction 
between the two species. In Fig. ^c) we show the re- 
sponse of the PSF phase. Due to the pairing order, both 
species move instantaneously and the time evolution of 
the center of mass is identical. The impulse on species 1 is 
transformed into a collective motion of both species and 
there is no relative motion. In the CFSF phase shown 
in Fig. m^d) , we see that the oscillatory motion of species 
1 is perfectly matched by an opposite or counter-flowing 
motion of species 2. We note that the sinusoidal oscilla- 
tion of the center of mass is not in contradiction with the 
constant current flow shown in Fig. [T] Averaging over 
lattice sites as in Eq. ([3]) erases this information. 

In Fig. [3] we plot the time evolution of the center of 
mass of both species, Xt{t), and the relative center of 
mass, Xr{t), after a brief displacement as in Fig. [2]for sev- 
eral interaction strengths. We increase the intra-species 
interaction U from 2 J to 20 J and U12 from 1.2 J and 12 J, 
with the fixed ratio U12/U = 0.6; we set Q/U = 10"^. 
The system undergoes a phase transition from SF to 
CFSF around U/J ^ 8. In the SF state, the brief dis- 
placement of species 1 leads to damped oscillations in 
both the relative and total center of mass. The maximum 



amplitude of the oscillation increases with U/J. The 
damping rate increases in this regime. As U/J increases 
beyond the critical point, the total center of mass mo- 
tion is suddenly suppressed to almost zero. On the other 
hand, the amplitude of the relative center of mass keeps 
increasing through the phase transition. The damped os- 
cillation of the relative center of mass suddenly changes 
to nearly undamped above the transition. 

We have also studied the dependence of the oscillation 
period on the interaction strength in the CFSF state. As 
noted previously, the density flow occurs at a constant 
velocity in this regime. The velocity can be estimated 
with two methods: the Luttinger liquid theory and an 
effective spin-1/2 model. Luttinger liquid theory predicts 
phonon-like low energy excitations with a velocity 
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AU12/J -BU12/U, 



(4) 



where the non- universal coefficients A and B are assumed 
to be constants. For a plateau of length L, the period of 
oscillation is then T — 2L/va- For the analysis of our 
simulations, we find it convenient to define the relation- 
ship based on Eq. (|4]), 
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l~ AU12/J - BU12/U. 



(5) 



In the limit of large U12/ J and U/ J, we can derive an- 
other relationship for the period. Assuming the particle- 
hole pairs are hard-core bosons, Eq. ([1]) can be mapped 
to an effective spin-1/2 model [8|. Using linear spin- wave 
theory [TtI, we obtain the phonon velocity va as 



VA 



Vi - C/12/C/, 



(6) 
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and 

(^)' = 4(JM2r(l-f/W^). (7) 

Unlike for the Luttinger theory, this relationship does 
not have free parameters. Figure ^c) shows the quan- 
tity {hLf/{dJTf as a function of U12/J for fixed U12/U, 
where L and T are obtained from the simulations. Over- 
all, the quantity {hL)'^ / (dJT)'^ shows a linear trend con- 
sistent with the relationship in Eq. ([5]). The value is in 
good agreement with Eq. ([7]) especially for U12/J > 8.4. 

In Fig. |3]we consider the time evolution after a sudden 
displacement for the SF, PSF and CFSF initial states. 
The displacement is the same for both species. For the 
SF and PSF initial states, the centers of mass of both 
species oscillate in phase around the new minimum of 
the trap. The origin of the in-phase motion is different 
for these two initial states. For the SF state, the species 
respond independently to the same displacement. For the 
PSF state, the species respond as pairs. The oscillation 
frequencies for the SF and PSF states differ. This is not 
only because the states are confined by different trap- 
ping potentials, but also because of the larger effective 
mass of the pairs in the PSF state. The oscillations ex- 
hibit damping that can be attributed to strong quantum 
fluctuations in ID systems [H [H El. For the CFSF 
state, the motion is overdamped and the center of mass 
remains at the original equilibrium position. The small 
oscillations that are visible in the time evolution of xi{t) 
and X2{t) are due to oscillations of the superfluid at the 
edges of the atomic cloud. For this type of displacement, 
the impact on both species is the same and it only affects 
the total density. Because the motion of the total density 
is suppressed for the CFSF state, the displacement can 
not induce a response. 

In Fig. [5] we consider a sudden displacement of the two 
traps in opposite directions. For the SF state, the two 
species oscillate around their respective new trap minima. 
The PSF state has no response, because the pairing order 



resists the force imparted by the trap displacement. The 
most dramatic response is in the CFSF state. Because 
the displaced traps act as an effective linear potential, 
2il{Di + D2)j, on the relative density, the two species 
are driven to move apart, until they reach the edge of the 
CFSF plateau. [3c)-(e) show density profiles at different 
times during the time-evolution. 

By increasing the displacement amplitude, we can po- 
tentially break the (anti-) pairing order. We can break the 
PSF pairing by applying a large opposite displacement, 
while the CFSF pairing can be broken by a large identical 
displacement. The displacement amplitude then becomes 
an indicator of the binding energy of the pairs in the PSF 
state and the anti-pairs in the CFSF state. 

In conclusion, we have studied the dipole oscillation 
of ID two-component Bose mixtures in an optical lattice 
and a weak confining harmonic potential. The oscilla- 
tion is induced by displacing the harmonic confinement 
either briefly or suddenly. We have shown that the re- 
sponse of the system shows much richer features than 
its single-component counterpart. In the two-component 
system, there are three long-range orders: the superfluid, 
paired superfluid and counter-flow superfluid orders. For 
the PSF and CFSF states, the suppression of individual- 
ity leads to distinct dipole oscillations dependent on the 
character of the perturbation. We have shown that this is 
the consequence of pairing and anti-pairing ordering for 
the two states. For the CFSF state, which forms plateaus 
of constant total density, the dipole oscillations resemble 
the oscillation of a conflned homogeneous system with 
free propagation between boundaries. 
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